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We recently reported that hairpin (or stem-loop) priming is better-suited than polyA tailing to generate cDNA for plant 
microRNA qPCR. One major limitation of this method is the need to perform individual cDNA synthesis reactions for the 
reference gene and test miRNAs. Here, we report a novel fusion primer that allows multiplexed hairpin cDNA synthesis of 
the most-commonly used reference gene, nucleolar small RNA U6, together with test miRNAs. We also propose the use 
of miR1515 as a house keeping control for tropical legumes. We show that multiplexed cDNA synthesis does not result in 
loss of sensitivity and reduces the amount of RNA required for miRNA gene expression assays. 



microRNAs (miRNAs) are short non-coding RNAs that regulate 
gene expression at the post-transcriptional level in both plants 
and animals (reviewed in refs. 1—4). In plants, their biogenesis 
starts with pri-miRNA transcripts which are capped and gen- 
erally poly-adenylated and then cleaved by DICER-LIKE 1 to 
give rise to short duplex RNA molecules consisting of the mature 
miRNA (typically 21-22 nt in length) and miRNA* molecules. 
Both miRNA and miRNA* are then 2' O-methylated at the 
3' end. Mature miRNAs are incorporated into RNA-induced 
silencing complexes and direct inhibition of translation or cleav- 
age of complementary mRNAs (reviewed in refs. 5—7). Reliable 
quantification of individual mature miRNA levels is necessary 
to understand their biological roles. Hairpin/stem-loop prim- 
ing** and PolyA tailing' are the two most commonly used cDNA 
synthesis methods for plant miRNA qPCR. We recently showed 
that hairpin priming method is better suited than polyA tailing 
method to synthesize cDNA for plant miRNA qPCR.'° One dis- 
advantage of hairpin priming method is the use of specific RT 
primers for the cDNA synthesis and therefore individual cDNA 
synthesis reactions for each miRNA. This not only requires more 
input RNA, but also results in possible differences in cDNA 
synthesis efficiencies between test miRNAs and reference genes. 
It is crucial that the reference genes and test miRNAs undergo 
identical reaction conditions during cDNA synthesis and qPCR. 
This ensures equal efficiency of the reaction between the refer- 
ence gene and the test miRNA(s) leading to accurate and reliable 
measurements. 



In general, other small endogenous noncoding RNAs such as 
U6, U24 or U26 are used to normalize the expression of both 
plant and animal miRNAs." '^ In Northern hybridization and 
polyA technique, the small nucleolar RNA U6 has been success- 
fully used as a normalization control for a wide range of experi- 
ments (e.g., see refs. 13 and 14). However, no suitable method 
exists for U6 cDNA synthesis coupled with hairpin priming. To 
resolve this issue, we designed cDNA synthesis and qPCR prim- 
ers for U6 from a highly conserved region of the gene (Fig. SI A). 
For U6 cDNA synthesis primer design, we added the "universal 
reverse" primer used in hairpin cDNA qPCR to the 5' end of 
U6-specific reverse primer sequence (Fig. SIB). When this modi- 
fied U6 reverse primer is used for cDNA synthesis, U6 expression 
can thereafter be assayed using hairpin cDNA qPCR technique, 
using U6-specific forward primer and universal reverse primer. 
"Hairpin-primed" cDNAs were generated using U6-specific 
cDNA synthesis primer from RNA samples obtained from dif- 
ferent plant species. RNA preparations were obtained from rice 
(leaf and seed), wheat (leaf), Arabidopsis (seedlings), medicago 
(leaf), corn (leaf meristematic tissue), grape (root), soybean (root) 
and prairie cord grass (rhizomes). All cDNA synthesis reactions 
were performed on an ABI thermocycler (GeneAmp 9700). All 
qPCR assays were performed using a MX3000P thermocycler 
(Stratagene/Agilent Technologies) and SYBR Advantage qPCR 
premix (Product* 639676, Clontech). Three different dilutions 
(1/50, 1/500, 1/5,000) of cDNA from each plant species were 
examined for U6 expression by qPCR (Fig. 1). Suitability of 
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Figure 1. Validation of the novel U6 fusion primer for cDNA synthesis on multiple plant species; (k) Arabidopsis, (B) soybean, (C) medicago, (D) grape 
root, (E) corn, (F) wheat, (G) rice leaf, (H) rice seed and (I) prairie cordgrass. Amplification plots of U6 in different plant species were obtained using 
primers designed in this study to adapt the use of U6 as normalization control in hairpin cDNA qPCR assays for miRNAs. Linearity (Rsq) values obtained 
by examining expression in different dilutions (circle: 1/50, triangle: 1/500 and diamond: 1/5,000) of cDNA are indicated in each panel. 



these primers to synthesize cDNA and assay U6 expression were 
determined by examining the linearity of amplification, quality 
of the amplification curves and dissociation curves. We tested a 
range of monocotyledonous and dicotyledonous (both legumes 
and non-legumes) species. The primers we designed were able 
to successfully detect U6 expression with reliable amplification 
efficiencies and linearity in different plant species (Fig. 1). 

One of the legume-specific miRNAs, miR1515 is uni- 
formly expressed in a variety of tissues in soybean" as well as 
other tropical legumes such as common bean and peanut."" To 
utilize miR1515 as a reference gene for miRNA qPCR in soy- 
bean, we designed hairpin primers for this miRNA. Expression 
of miR1515 was assayed in different soybean tissues and in the 
roots during nodule development (Fig. 2). For expression analysis 
in different soybean organs, samples were collected from root, 
stem and young leaves of 12-d-old plants, 14 d post inoculation 



(dpi) mature nodules, flowers and young pods. For time course 
experiment with soybean, root tissues were harvested at 0, 1, 3, 
6, 12 and 24 h and 3, 7 and 14 dpi from mock- (nitrogen free 
plant nutrient solution) and B. japonicum-'moculaxed plants. All 
root tissues harvested were first washed thoroughly, blot dried, 
immediately frozen in liquid N, and stored at -80°C until RNA 
isolation. RTqPCR analysis showed the uniform expression of 
miR1515 in the different tissues (Fig. 2A) and along the course 
of nodulation (Fig. 2B), as previously identified using Northern 
analysis." We propose the use of miR1515 as reference gene for 
miRNA expression assays in soybean and possibly other tropi- 
cal legumes (e.g., common bean, cowpea, peanut). The use of a 
miRNA as reference gene for miRNA expression assays has some 
advantages over using protein coding genes such as actin, ubiqui- 
tin or EFl. Since miRNA biogenesis pathways are distinct from 
mRNAs, alterations in biosynthesis machinery (e.g., due to virus 



e24918-2 



Plant Signaling & Behavior 



Volume 8 Issue 8 



o 
o 



2 
.o 



y) 

I— 

Q. 

X 
(D 
0 

> 

0) 



B 



.18 

0 o 
0 c 

> o 

i5 CD 
0 '^-^ 



infection) might introduce artificial abundance differ- 
ence between miRNAs and mRNA reference genes used 
as normahzation controls. Using a miRNA as reference 
gene can help overcome such issues. 

Finally, we examined if cDNA synthesis for U6 and 
miR1515 could be multiplexed along with test miRNAs 
without loss of efficiency. Hairpin primers for up to four 
different miRNAs (50 nM each final concentration) and 
the reference gene (U6 or miR1515) were added to a sin- 
gle cDNA synthesis reaction using soybean root RNA as 
template. Fold difference in miRNA expression between 
samples or between miRNAs and respective reference 
genes were calculated using the ddCt method.'^ We com- 
pared Ct values between qPCR assays performed using 
multiplexed cDNA or independently synthesized cDNA 
(Fig. 3). There was no significant difference in efficiency 
or Ct values when miRNA abundance was assayed in 
single miRNA cDNA synthesis reactions compared with 
multiplexed cDNA synthesis reactions (ANCOVA; p > 
0.05). Therefore, multiplexing during cDNA synthe- 
sis would not only enable the efficient use of minimal 
amounts of RNA for qPCR, but also enhance the accu- 
racy of the assays by ensuring equal cDNA synthesis effi- 
ciency between test and reference genes. 

There are other situations where multiplexing in hair- 
pin cDNA qPCR might be handy. The ability of hairpin 
cDNA qPCR to distinguish family members or mature 
sequence variants can sometimes be a limitation, e.g., 
when attempting to assay miRNA abundance of mul- 
tiple/all family members. Nevertheless, this could be 
addressed by multiplexing cDNA synthesis for multiple 
family members /variants in a single reaction using spe- 
cific hairpin primers for each variant/family member. 
This limitation can also be addressed by the use of a 
recently developed method (miQPCR, reviewed in ref 
18) where an RNA adaptor is ligated to the 3' end of 
mature miRNAs and the ligated molecules are reverse 
transcribed using a primer complementary to the adap- 
tor sequence. However, this approach might as well suffer 
from the inability to distinguish methylated vs. non-methylated 
RNA molecules, especially in RNA preparations from plants. 

In this study, we have developed tools to enhance hairpin 
cDNA qPCR method to assay plant miRNA gene expression. 
First, we developed a cDNA synthesis primer for U6 that allows 
multiplexed cDNA synthesis and validated it using multiple 
plant species. Next, we developed miR1515 as a reference gene 
to normalize miRNA expression in soybean and possibly other 
tropical legumes. Finally, we successfully multiplexed cDNA 
synthesis of multiple miRNAs with reference genes without loss 
of efficiency. 
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Figure 2. Expression of miR1515 in different soybean tissues (A) and along the 
course of nodulation (B) in soybean assayed by hairpin cDNA qPCR. (A) Relative 
expression compared with roots were calculated using average Ct values from 
four independent replicates. (B) Relative expression at each time point post 6. 
japonicum inoculation compared with mock inoculated roots. Data shown are 
average of four independent biological replicates. Error bars represent SD. 
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Figure 3. Multiplexing during hairpin cDNA synthesis did not affect assay efficiency. Ct values from hairpin cDNA qPCR assays to quantify (A) miR1515 
(B) miR164 and (C) miR166 plotted using different dilutions of cDNAs obtained through independent (filled circles), duplex (open squares) or multiplex 
(filled triangles) synthesis reactions. Error bars indicate SD. In all cases, the slope of the curves or Ct values (at each dilution) did not differ significantly 
between different cDNAs. 
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